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Abstract 

An analytic ground state is proposed for the unbiased spin-boson Hamiltonian, 
which is non-Gaussian and beyond the Silbey-Harris ground state with lower ground 
state energy. The infrared catastrophe in Ohmic and sub-Ohmic bosonic bath plays 
an important role in determining the degeneracy of the ground state. We show that 
the infrared divergence associated with the displacement of the nonadiabatic modes in 
bath may be removed from the proposed ground state for the coupling a < ac- Then ac 
is the quantum critical point of a transition from non-degenerate to degenerate ground 
state and our calculated ac agrees with previous numerical results. 
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I Introduction 



Quantum impurity systems with competing interactions constitute a field of wide in- 
terest in the quantum physics. In recent years, the quantum two-level system coupled 
to dissipative bosonic environment (spin-boson model, SBM) attracts much attention 
in this field because it may be one of the simplest but nontrivial quantum impurity 
system for studying the physics of competing interactions. The Hamiltonian of SBM 
reads (we set h = 1) 

H = -^^(^x + J2ujkbibk + ^Y.9k{bl + bk)a,, (1) 

k k 

where b^ {bk) is the creation (annihilation) operator of environmental bosonic mode 
with frequency Uk-, and o"^ are Pauli matrices to describe the two- level system. The 
competing interactions in SBM are between the quantum tunneling A and the dissipa- 
tive coupling Qk to the environment. The effect of the environment is characterized by 
a spectral density JioS) = J2k 9k^{^ ~ '-^k) = 2au'^ul~'^9{ijJc — (jj) with the dimensionless 
coupling strength a and the hard upper cutoff at Uc- The index s accounts for various 
physical situations[Tl |2]: the Ohmic s = 1, sub-Ohmic s < 1 and super-Ohmic s > 1 
baths. 

The quantum critical point (QCP) and the quantum phase transition (QPT) are 
related to the ground state transition, which is usually triggered by competing inter- 
actions. As for SBM, the interesting phase transition is related to the transition of 
degeneracy of the ground state, that is, it is a transition between the non-degenerate 
and degenerate ground state [H [21 |3l HJ |5] . The main theoretical interest of the QCP 
in SBM is to understand how the competing interactions inffuences the degeneracy of 
the ground state. Since the Hamiltonian (1) is invariant under cr^ — )■ — (together 
with bk,bl —7- —bk, —bl) and one must have {o'z)g = ((...)g nieans the ground state 
average). However, for the Ohmic bath s = 1 it is well knownjll [2] that a Kosterlitz- 
Thouless quantum transition separates a degenerate ground state at a > etc from a 
non-degenerate one at a < (c^c = 1 in the scaling limit A ^ Uc)- 

The ground state of SBM Hamiltonian (1) was studied by many authors using 
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various analytic and numerical methods. Silbey and Harris (SH)[6] proposed an varia- 
tional ground state and predicted the QCP etc = 1 for s = 1. The SH ground state was 
used by Kehrein and Mielke[7j for sub-Ohmic {s < 1) bath to calculate the QCP 
In last ten years, various numerical techniques were used for calculation of the QCP 
in the SBM, such as the numerical renormalization group (NRG)[3l HJ |5], the quan- 
tum Monte Carlo (QMC)|8j, the method of sparse polynomial space representation [9], 
the extended coherent state approach[TO], and the variational matrix product state 
approach Besides, recently an extension of the Silbey-Harris ground state was 
proposed by Zhao et a/. [12] and Chin et a/. [13] to study the QPT in the s = 1/2 
sub-Ohmic SBM. 

In this work, we propose an analytic ground state wavefunction for the SBM, which 
is non-Gaussian for the bath modes and is an extension of the work of Zhao et al.[12\ 
and Chin et al.[13\. The QPT is usually not a weak coupling problem and people be- 
lieve that the numerical techniques may be more powerful than approximate analytic 
methods for strong coupling problem. Then, why do we still try to find an approxi- 
mate analytical solution? Generally speaking, our purpose is to see and understand 
the physics more clearly and straightforwardly. In particular, here our purpose is to 
understand the role played by the infrared divergence in the SBM Hamiltonian (1). 

The QPT in quantum impurity systems may be related to the infrared catastrophe 
in baths. P. W. Anderson [14J was the first to point out this relation for the Anderson 
model and Kondo model in fermionic bath. Our question is: What is the role played 
by the infrared catastrophe in the quantum phase transition in bosonic bath of SBM? 

II The ground state 

If A = 0, Hamiltonian (1) is solvable and we have degenerate ground state 

l^ta)) = expl-Y^gM - bk)aj2uk]\ t (i))|{04), (2) 

k 

where | t (i)) is the eigenstate of cr^: cr^l '\ (J,)) = +(— )| t (i)) \ {^k}) is the vacuum 
state of the bath. Then, for finite A it is naturally to use a superposed ground state 
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to remove the degeneracy. But it is well known [H [21 [6l [7] that there exists an infrared 
divergence in the overlap between the degenerate states: (V'tlV'j,) = exp[— J2k dlf^^l] = 
for s < 1. Silbey and Harris proposed a modified superposed ground state [6] 

\Gsh) = exp[- J2 9kibl - h)aj2{uk + Vom^'^'il t) + I i))|{0.}), (3) 

k 

with finite renormalized overlap rjo = e^p[—J2k olf^i^k + ^o^)^] where the infrared 
divergence has been removed. The ground state energy is 

E^"" = -r]oA/2 -Y.9l{^k + 2r/oA)/4(a;fc + r/oA)l (4) 

k 

For the SH ground state at the scaling limit A ^ Wc, 770 = (eA/wc)!^ for s = 1 and 
thus the QCP is at af^ = 1 where r/o = 0. For sub-ohmic bath s < 1 one can calculate 
the QCP by condition: r^o = at a — )■ Q;f^[Tl [15], and some results are listed in the 
second column of table 1. 

Zhao et al\12\ and Chin et a/.[13j proposed an extension of the Silbey-Harris ground 
state to study the QPT in the s = 1/2 sub-Ohmic spin-boson model, with degenerate 
ground state when zero-biased and a > (superscript "D" means degenerate) 

= exp(-5±)(n±| t) + v±\ i))|{0,}), (5) 

^± = E ^(bl - hMkcrM^ - (6) 
k ^'^fc 

where u+ = = 2-i/Vl + M, u_ = v+ = 2-i/Vl - M, = Uk/{ujk + W), 
W = r^A/Vl -M2, and 



T] = exp 



T.9iek/2ul 



(7) 



M = J2glM^-^k)V{ookW). (8) 

k 

Zhao et al.[12\ and Chin et al.^3\ let 0^ = M to be a constant in Eqs.(6) and (8), thus 
|\E'+) and are degenerate with degenerate ground state energy 

= -W/2 - E&(2 - ^k)/^co, + Y.9lM\l - ikf/^^k. (9) 

k k 



Besides, the average {^±\<Jz\'^±) = may be finite. They proposed that the QCP 
is at a = where a nonzero M leads to lower ground state energy (Note that when 
a < , M = and |^+) = |^_) = \Gsh) )• Some values for different baths are 
listed in the third column of table 1. But, as mentioned above, since the Hamiltonian 
(1) is invariant under az — > — (together with bk,bl — )■ —bk,—bl) we should have 
{(Tz)g = 0. 

Ill The infrared catastrophe 

The wavefunction of every bath mode in \Gsh) or |\E'±) is a Gaussian function, thus 
these ground states are in the Gaussian approximation. Following the proposal of Shore 
and Sander[T6] we propose the following superposed ground state for the SBM, which 
is beyond the Gaussian approximation and takes into account the effect of quantum 
fluctuations, 

\G)=A{\^^) + \^.)), (10) 

where A is a normalization factor. Then, it is easy to check that {G\crz\G) = 0. But 
if one choose 0^ = M in Eqs.(6) and (8), as was pointed out by Chin et a/.[13j, there 
is an infrared divergence of the occupation number of the nonadiabatic (NA) modes. 
We show that this divergence leads to the orthogonality catastrophe between and 

P = = m}\ exp f- E ^(1 - - b,)M] |{0 J) 

-xp (- E - ^ exp [-cM^W^- I ^ 

for Ohmic {s = 1) and sub-Ohmic {s < 1) baths as the integration in the exponential is 
infrared divergent. This is similar to the infrared catastrophe in Fermi sea interacting 
with a quantum impurity|14]. Because of the orthogonality catastrophe the ground 
states, \Gd) = |^+), \Gd) = |^-) or \G) (Eq.(lO)) are degenerate with ground state 
energy (9). 



(ir 
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The way to avoid the infrared catastrophe is similar to the proposal of Anderson 
that is, quantum fluctuation of the NA modes leads to a /c-dependent 0^ in Eqs.(6) 
and (8) removing the infrared divergence. Then the ground state energy Eg of the 
superposed ground state (10) is 



Eg = {Eo + pU)/{l + pVr^M2), 
where Eq = p = (^+|^_) and pU = Here 

Eo = -W/2 - ^ glW^ - a)/4u;fc + Y, 

k 

\ k f 

-riA [cosh(M) - 1 - M(sinh(M) - M)] /2, 

andy = Efc^7M(l-a)V4c^fc- 

The variational function 0^ can be determined by 



(12) 



(13) 



(14) 



dEg dE^m ^ dE^dp_ ^ dEg dY 



DM 



dp d(f)k dY d(j)k 



for every mode k and the result is 



0fc = ruk/{ujk + pS), 



(15) 



(16) 



where 5 = 2{EoVl - - U)/[{1 - p)(l + pVl - M'^)] and r is the variational pa- 
rameter. In this way, the overlapping integral is 

;i - 4fc)'0fc j = exp l^-ar^W^ [ j—^^ '^'^ 

which is finite as long as s > 0. 



p = exp 



9l 



[u + wy{uj + psy 



,(17) 



For s = 1 the result of variational calculation is shown in Fig.l. When a goes to 
1, the variational parameter r tends to 1 and the overlapping p decreases to zero as 
follows: 



6_ 
W 



exp 



1 — ar^ 
6 



ln(l + W/u,) + 



2 + W/uc 
1 + W/oOc 



18) 



that is, p — )■ when a — 1 — 0^ (0+ is a positive infinitesimal) since r — 1. This is to 
say that for s = 1 the ground state becomes doubly degenerate when a — )■ = 1. 

etc for sub-Ohmic bath (s < 1) can be calculated in the similar way, that is, the QCP 
where the ground state changes from non-degenerate {a < etc) to doubly degenerate 
{a > etc). Our results for some s values are shown in Table 1. For comparison, the 
numerical results by NRG[3J, by QMC[8], by the method of sparse polynomial space 
representation |9], and by the extended coherent state approach[10] are also shown. One 
can see that our result compare well with these numerical results. 

Fig. 2 shows the difference between our calculation of the ground state energy and 
that of Zhao et a/. [12] and Chin et a/. [13], SEg = Eg — E^. The lower ground state 
energy indicates that the ansatz of this work is a better one for the real ground state. 

We note that when s > 1 (super-Ohmic bath) the overlapping p in Eq.(17) has 
always a finite solution. This is to say that the ground state of the SBM with super- 
Ohmic bath is always non-degenerate and there is no QPT. 

IV Conclusion 

We propose an analytic ground state wavefunction for the unbiased spin-boson Hamil- 
tonian, which is a superposition of the two degenerate state and is non-Gaussian for 
the bosonic bath modes. The infrared catastrophe in Ohmic and sub-Ohmic bosonic 
bath plays an important role in determining the degeneracy of the ground state and 
we show that the infrared divergence associated with the displacement of the nonadia- 
batic modes in bath may be removed from the proposed ground state for the coupling 
a < Oc- The QCP ac is determined by the transition from non-degenerate to degener- 
ate ground state. Our ground state energy is lower than previous authors'results. The 
calculation of ac agrees well with previous numerical results. 
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Figure 1: The variation parameter r and the overlapping p as functions of a for Ohmic 
bath s = 1. 
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Figure 2: Eg — E^ is the difference between ground state energies calculated by Eq.(12) 
and Eq.(9). A/wc = 0.1. See text for details. 
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Table 1 QCP of different bath type s. 
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etc [3J 


a,[8\ 
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1/4 


0.08554 


0.02413 


0.02744 


0.0264 


0.0254 


0.0259 


0.0256 


1/2 


0.1768 


0.08555 


0.1084 


0.1065 


0.0983 


0.0977 


0.0 


820 
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0.3076 


0.3168 


0.2951 


0.2953 


0.3205 
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0.5121 
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1 
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